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We study a class of exactly solvable models for strongly 
correlated electrons, defined on a set of A*' cells, and with in- 
finite on-site repulsion on part of the sites of each cell. For 
2N or more electrons the exact ground state is known. We 
construct a tractable (2A'' — l)-particle state which becomes 
asymptotically degenerate to the ground state in the thermo- 
dynamic limit for one special D-dimensional model. For other 
models, that state may be used to calculate variational upper 
bounds on the ground-state energy. For a Hubbard chain with 
three sites per unit cell, the analytical bound is compared to 
numerically exact results. 
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I. INTRODUCTION 

The theoretical description of electronic correlations in 
solids is still a challenge. Even one of the simplest models 
for correlated electrons, the Hubbard model,EJ exhibits 
enormous mathematical complexity. Its analysis is far 
from being complete. 

Recently, a new class of models for tight-binding elec- 
trons with infinitely strong on-site repulsion has been dis- 
covered. For particle numbers of two or more particles 
per unitjjell, one can construct the ground states of these 
models.H These ground states are spin-singlets .and have 
the structure of resonating-valence-bond stateaa for cer- 
tain cases. 

The class of models which are solvable -biz. this 
method has been generalized by-several authors .0~Q One- 
dimensional models of this classO allowad for the calcula- 
tion of equal-time correlation functions^ using a transfer- 
matrix technique. All equilibrium correlations studied so 
far decay exponentially. Dynamical correlations obtained 
numerically for the one-dimensional Hubbard chain in- 
dicate dispersionless excitation spectra for an electron 
density of two per cell.El The propagation of one hole, 
however, shows dispersive delocalized features. 

The regime of particle numbers below two per unit cell 
has not been accessible by analytical methods so far. The 
main result of this paper is an analytical upper bound 
on the ground state energy of systems containing 2N — 1 
particles {N denotes the number of cells). The trial state 
and the corresponding upper bound on the ground-state 
energy are discussed for two examples. The first example 
is a hypercubic Hubbard model as introduced in Ref. |^. 
It is shown that the trial state becomes an asymptotically 



exact eigenstate of this particular (2A^ — 1) -particle sys- 
tem in the thermodynamic limit. For the second system, 
a linear Hubbard chain, the upper bound on the ground- 
state energy may be calculated analjrtically as well. A 
comparison to exact numerical resulta3 is presented. 

II. THE CLASS OF SOLVABLE MODEL 
HAMILTONIANS 

In this section we recall the_most general description 
of the class of solvable models.El This class contains Hub- 
bard, Anderson and Emery models in the limit of in- 
finitely strong interaction (U = oo). The solution pre- 
sented below requires certain lattice structures. (In the 
following, the term graph is used, because lattice implies 
translational invariance which needs not be assumed.) 

The models are defined on graphs, where each of the 
N vertices contains a cell of sites. (The cells need not 
be identical either.) The cells are labeled by an index 
i = 1, ... TV. Let the set of all sites within the cell i be 
called Afi. (It is mentioned in passing that A/i is not a 
unit cell if the graph is translationally invariant, because 
neighboring cells may share some sites.) The sites (or 
electronic orbitals) within a cell i are labeled by an index 
a. Some of those orbitals (not necessarily all) may carry 
an infinitely strong on-site repulsion, e. g. a cell may con- 
tain "d-sites" which may be empty or occupied by one 
particle only (with spin up or down) and "p-sites" which 
can be occupied by up to two particles. The set of d-sites 
is called Ui. The repulsion is incorporated into the model 
by a projection operator 

p = n n (1) 

which strictly excludes double occupancy on all sites 
«,Q! G Ui, where n^^^^^ = cl^ ^^Cj^^ g.. The Fermi oper- 
ators c\ ^ ^ i'^i a a) denote creation (annihilation) opera- 
tors and are defined as usual. The connection between 
cells is established by common sites. We denote the set 
of common sites by the intersection dj — J^fiCi N'j. Ex- 
amples for some models are given in Ref. §,§,0. Figure |l| 
illustrates the two cases discussed below. 

In order to construct the ground state of the models 
defined above consider the linear combination of Fermi 
operators 

*U --^ E Ko^Ac.^ ■ (2) 
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They obey the following algebra: 

K.'*L']+-[*..'*jv]+ = (3) 
l-^. vfit 1 ^ J 

f E WaJ^ for 

I E '5(»a),(i/3)^U,^^,/3.<T Otherwise 

(The Kronecker symbol (^fj^) q-^) assumes the value one if 
the indices ia and j(3 label the same site and equals zero 
otherwise.) The coefficients Aj q, g. define the parameters 
of the Hamiltonian which will be constructed below. 

The \['-onerators allow to write down a positive 
semidefinitdHI Hamiltonian in a very concise way: 

H = J2Ha with H, = Y,^,^,P^l^. (4) 

cr i 

However, the physical meaning of this Hamiltonian be- 
comes more transparent if it is rewritten in terms of the 
Fermion c-operators by using the relation 

"fLP = + E >^^.aA»,an^.o.,-aP (5) 

aeu, 

and the definition of the operators as defined in Eq. (||). 
The Hamiltonian (^) is transformed into the following 
form: 

^cr = - P \ E K.a,crK,0.a (H,a,a^i,P,a- (6) 

+ E IWPK„.^-i) 

+ E \Kc.J'n,^».-a}P- 

The first term in describes particle hopping within 
a cell where each site (orbital) is connected to all oth- 
ers. (Particle transfer between interacting orbitals and 
orbitals without interaction is usally called a hybridiza- 
tion. In that case, the model would describe e. g. an 
Anderson model.) The remaining terms contain trivial 
constants and sums over occupation numbers which add 
up to the total particle-number operator for the case of 
the I?-dimensional Hubbard model as discussed below. If 
the sets Ui and J\fi are not identical or if not every site 
within a cell is connected to a neighbor cell, the sums 
over occupation number operators may add up to an ad- 
ditional field which shifts the energy of certain sites. Such 
a model then contains more than one parameter and the 
exact ground state as derived below is only valid on a 
parameter surface defined by the quantities A. |Ih the 
following, all A-coefficients are assumed to be realty. 



III. THE GROUND STATE FOR 2N OR MORE 

PARTICLES 

The models as defined above allow to write down an 
exact ground state if the system contains two or more 
particles per cell. The ground state is a Gutzwiller pro- 
jected Slater determinant: 

|a>o) :=Pn*lT*lil^)- (^) 

i 

The expression i?|$o) always contains a factor (P'^l ^) ■ 

Because of the equation {P'^l ^) — 0, the state Eq. (0) 
is an eigenstate of H (^), the corresponding eigenvalue 
equals zero. r-. 

Further, H is positive semi-definitecl. Therefore the 
state (0) belongs to the ground-state manifold. In Ref. 
Tasaki provides a proof of the uniqueness of the ground 
state if the state \x) is equal to the vacuum state |0) (27V 
particles). Larger fillings can be discussed if \x) contains 
additional particles, for that case it immediately follows 
by construction that the state (^ is degenerate. 

If all sites carry an infinitely strong on-site repulsion 
and, e. g. the A-coefficients are set to unity, the state 
(0) exhibits_a structure often called resonating-valence- 
bond state.EJ In that case the ground state consists of a 
linear combination of local singlet bonds jj P^'j^^^'^^ = 

P Ea#/3 {4,a,fU,i + 4/3,T4a,i) ■ projcctor ensures 
that no two singlet bonds have a site in common. 

IV. UPPER BOUND ON THE GROUND-STATE 

ENERGY FOR 2N-1 PARTICLES 

Although the ground state is known for two or more 
particles per cell, the only analytical result that has been 
obtained for particle numbers Ne below the "magic" 
number — 2N so far is the trivial lower bound on 
the ground state energy Eo{Ne) > which is valid for 
any number of particles. [The Hamiltonian (^ is non- 
negative.] The present section provides an upper bound 
on the ground state energy of H for one additional hole, 
e. g. a particle number of iVe = 2A'^— 1 using a variational 
state. The quality of variational methods is not very well 
controlled, there exists no small parameter, only an in- 
equality for the energies. Therefore the method requires 
some intuition in order to "hit" the correct physics by 
guessing a good state. The first guess of a variational 
state, the state 

1^') ^=^*It n *It*lio), (8) 

j, i=^i 

is led by the idea to leave as much of the 2A^ -particle 
physics unchanged as possible. The state is constructed 
in a way that most terms of the Hamiltonian generate 
zero-states due to the relation [P^Xaaf' ~ 0- From 
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= and Hi\ipi) = *i.J$o) (|*o) denotes the 
ground state of the system with 2N particles of the pre- 
vious section) one obtains the diagonal matrix 



(9) 



If the norm of the states \Lpi) and |$o) were known, 
the expectation value of H could be evaluated. Unfor- 
tunately, this is not possible in general. However, the 
following orthogonal and normalized [see Eq. (||)] trial 
states 



H 



(10) 



allow us to calculate the desired expectation value. (The 
positive square root of H is well defined because H > 0.) 
Surprisingly, the Hamiltonian matrix elements {Xj\H\xi) 
of the {2N — 1) -particle system can be expressed in 
terms of occupation-number expectation-values of the 
2iV-particle system in the following way: 



^ aeJ\fi 



(11) 



J2 \.Q, 

aeu^ 



a,l3£{CijnUi} ' 



where 



($o|*o) 



(12) 



In the following section, the upper bound Eq. (|ll|) on the 
ground-state energy is discussed for two examples which 
allow us to calculate the expectation value (rij ^ ^). 



V. RESULTS 

A. D-dimensional Hubbard model 

The first example, where the upper bound on the 
ground-state energy of the (2A^— 1) -particle system may 
be evaluated is the Hubbard model as introduced in Ref. 
^, where all A's are set to unity and the spatial dimension 
is larger than one. For the three-dimensional case, the 
cells have the topology of octahedra with 2D ~ 6 sites. 
All sites carry interaction, therefore the sets Ui and Mi 
are identical. Further, every site in Mi is shared by ex- 
actly one neighboring cell, forming a D dimensional (hy- 
per)cubic lattice of cells. Periodic boundary conditions 
are assumed. The corresponding lattice is visualized in 
Fig. I^a. The Hamiltonian of this system reads 



/ J 1.(7 l.O 



(13) 



2D 



i.i7 a,/3— 1 



The last term in (|13|) contains a trivial constant 'iND and 
the total electron-number operator N^. (with eigenvalue 
Ng) which is a conserved quantity. 
With the Fourier transform 



\Xk) 



,ikj 



Xj) 



(14) 



and (fti^Q,_(j) = D^^ (the occupation probabilities for all 
sites are eawal for symmetry reasons) one obtains an up- 
per boundlla on the (2A^ — 1) -particle ground state energy 



D 

{Xk\H\xk) = 2(1 - -^){D + ^ cos fc,0 



D' 



(15) 



The minimum Ey{2N — 1) = imn {xk\H\xk) with re- 
spect to the fc-vector is obvious: For (fci, /c2, ■ • ■ , ko) = 
(tt, TT, . . . , tt) the right hand side of ([l5| ) becomes zero. 
However, this fc-vector has to be exdnded for this partic- 
ular model, because it can be showna (c. f. App. ^ that 
the state |<I>o) (0) vanishes if fc = (tt, . . . , tt) is allowed. 
That particular fc-vector can be excluded by the restric- 
tion to systems where one of the D system dimensions 
contains an odd number of cells. Then the fc-vector com- 
ponent corresponding to this particular direction cannot 
equal TT. 

This result indicates that the state \xk) becomes an 
asymptotically exact (27V — 1) -particle ground state of 
H for the thermodynamic limit. The upper bound on 
the ground state energy vanishes as 



Ey{2N ~ 1) = 0{N-1^) 



(16) 



and asymptotically approaches the lower bound i?< {2N— 
1) = 0. This result is even relevant for D = 2, where 
the asymptotically exact state describes one hole in a 
half filled RVB-background. The absence of_ferromag- 
netism does not contradict Nagaoka's theorenJl3 because 
the present lattice is not bipartite. (For a bipartite lat- 
tice, the state would be a ferromagnet, i. e. the total 
spin would assume its maximum value.) Unfortunately, 
no rigorous statements on the uniqueness of this ground 
state can be made. 

Further, it can be seen that for D > 2 the chemical 
potential of the system for particle numbers of Ne — 
2N — 1, . . . , ND remains constant. The non- interacting 
tight-binding analogue exhibits quite similar behavior: 
The band structure of this particular lattice consists of 
one cosine band with e(fc) = 2 — 2{D + '^^j^cos k^j^) and 
D — 1 flat bands at the upper band edge of the extended 
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band. However, the "bandwidth" of the interacting sys- 
tem is reduced by a factor (1 — 1/D), which may be in- 
terpreted as an increased mass. This result sheds new 
light on the conductivity in the system. It might be 
that the many-particle state describes a dispersive delo- 
calized hole which may contribute to conduction. Unfor- 
tunately, no results on hole densities have been obtained 
and the dispersive state is not exact except in the vicinity 
of fc — (tt, . . . , tt). 



B. Linear Hubbard chain 

As a second example, we study a linear chain with N 
cells and 4 sites per cell A/i (3 sites per unit cell). The 
topology of this chain is illustrated in Fig. |l|b; it may be 
interpreted as a chain of connected tetrahedra. Again, 
periodic boundary conditions and the case of infinite on- 
site repulsion on every site are assumed, and thus the sets 
Afi and Ui are identical. The cells i and i+1 are connected 
by one site (backbone site): the sites with indices (i,a = 
1) and {i + 1, f3 = A) are identical. All A's are set to unity. 
The Hamiltonian of this particular chain reads 



(17) 



i,cr cf3=l 



+ P (^8N - 27Ve) . 

In addition to the hopping terms, one obtains a local field 
which shifts the energy of the backbone sites {a = 1). 

For this model, the_npper boimd Eq. ( pl| ) reduces to 
the simple expressiorO 

{Xk\H\xk) =3 - (n,^„^i_|) -f 2(1 - (n^^^^i |))cosfc 

(18) 

The right hand side of Eq. (^|) assumes a minimum 
for k — TT (which is allowed for this model) . This leads to 
the following upper bound on the ground state energy: 



Ey{2N - 1) = 1 



(19) 



The occupation probability on the backbone sites may 
be calculated by a transfer-matrix technique for the ther- 
modynamic limit. This calculation is briefly outlined in 
appendix H and yields 



lim Ey(2N- 1) = -(3 ^ 



) w 1.361325 (20) 



In contrast to the previous example this result indicates a 
gap in the (2A^— 1) -particle spectrum. Equation ( pO| ) es- 
timates this gap from above. Therefore the upper bound 
is compared to numerically obtained exact results as de- 
scribed in Ref. pi 



The numerical results hajie-Jpeen obtained by the 
conjugate-gradient algorithniEj'tj The energies are listed 
in Tab. |. For chain lengths up to four unit cells, the 
ground state is calculated without applying symmetry 
restrictions. In order to check for ground-state degen- 
eracy, the determination of the lowest eigenstate is re- 
peated in the subspace orthogonal to the ground state. 
For an even number of unit cells, the energies of these 
orthogonal eigenstates do not coincide with the ground- 
state energy, therefore the ground states for even chain 
length are non-degenerate. (Chains with an odd number 
of unit cells are degenerate because of reflection symme- 
try which induces degeneracy with respect to ±fc.) 

As a second step, the same calculations have been per- 
formed assuming site-interchange symmetry as decribed 
in Ref. ^ and translational invariance with k = tt. The 
energies of the symmetrized states coincide with the pre- 
viously calculated non-degenerate ground-state energies, 
therefore the symmetry of the ground state is uniquely 
determined. 

The chain with six cells is numerically treated by appli- 
cation of site interchange symmetry only and the result 
for A'^ = 8 is obtained assuming translational invariance 
in addition to site-interchange symmetry. 

The numerically obtained exact ground-state energies 
Eo{2N — 1) are shown in Tab. |. The upper bound Eq. 
( pO| ) of the infinite system deviates from the numerical 
result for A^ = 8 by « 13%. This agreement is not very 
close. However, the symmetries of the variational state 
and the numerically exact state coincide. Better quan- 
titative agreement may be obtained by the expectation 
value of H with respect to the state 



(21) 



[The state \ipi) is defined by Eq. (^).], which may be eval- 
uated numerically. For a chain with A^ = 4 one obtains 
{ipk=TT\H\Lpk=-K) — 1.2555. This result deviates from the 
numerically exact result by « 4% only. The state \xk) 
is proportional to \/i7|</5fe), therefore a difference in the 
expectation values {(pk=7T\H\ipk=-K) and {xk=7v\H\xk=n) 
simply reflects the fact that |xfc=7r) (and hence \Lpk=Tr)) 
is not an eigenstate. The operator ^/H "amplifies" the 
deviation from the exact ground state. 

The single-particle band-structure of this particular 
chain is illustrated in Fig. || and consists of two extended 
bands, one broad band separated from a narrower one by 
a gap of the size of one energy unit. The highest energy 
band collapses to one dispersionless flat band which is 
again separated by a small gap. For 2A^ particles the 
non-interacting system would have a completely flUed 
lower valence band describing an insulator. It is obvi- 
ous that the interacting system shows different behavior 
here: In the framework of the one-particle language, one 
would expect a gap if the particle number is changed 
from A^e = 2A' to Ne = 2N + 1, however, for that case 
the chemical potential of the interacting system remains 
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constant. In contrast to the one-particle picture, the 
chemical potential of the interacting system changes if 
the particle number is changed from Ne — 2N — 1 to 
A^e = 2A^. This contradiction illustrates the breakdown 
of one-particle physics in this particular model. 



VI. SUMMARY 

For a class of models describing strongly interacting 
Fermions it is possible to write down the exact ground 
state if the particle number equals 2 per unit cell or 
more. The particle number regime below this threshhold 
has not been accessible by rigorous methods so far. We 
present a tractable variational state for systems with one 
additional hole. The corresponding expectation value of 
the Hamiltonian approaches a lower bound for the case of 
a particular D-dimensional Hubbard model in the ther- 
modynamic limit. We conclude that the presented varia- 
tional state approaches the ground state or an eigenstate 
which is degenerate to the ground state if the ground 
state is not unique. 

For a second model, a special Hubbard chain, the an- 
alytical expectation value deviates significantly from the 
numerically calculated ground state energy. Although it 
is shown that the numerical ground state and the varia- 
tional state exhibit the same symmetries, it remains un- 
clear, how well the ground-state physics of the chain is 
described by our variational ansatz. 
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APPENDIX A: 

This appendix provides a proof that the state |$o) [Eq. 
(^] of the 2N- particle D- dimensional system vanishes, 
if the system extends an even number of unit cells in each 
direction. 

The state and the Hamiltonian are expressed in terms 
of creation (annihilation) operators (^i o-) referring 
to non-orthogonal one-particle states. However, they 
must be linear independent, otherwise the correspond- 
ing Slater determinant vanishes. 

There is no problem for models with non-equivalent 
sites, however, the _D-dimensional Hubbard model re- 
quires some care. Consider the Slater determinant 



The norm of 1$) can be expressed in terms of a determi- 
nant 



i*)=n*iT*iiio)- 



(Al) 



V($j¥) = det M 
The matrix M has the elements 



'<ij> J 



(A2) 



(A3) 



where 5,c.ij> equals one if i and j are nearest neighbors 
and zero otherwise. The structure of M is similar to 
a D-dimensional hypercubic tight-binding problem with 
periodic boundary conditions. 

The determinant is equal to the product over all eigen- 
values of the above tridiagonal matrix. The eigen- 
values are given by the familiar cosine bands 2{D + 
'Yl!^=i coskjj). If all system length contain an even num- 
ber of unit cells, one of the eigenvalues vanishes, be- 
cause one of all possible k- vectors equals (tt, tt, tt, . . . , tt). 
Therefore, the state |$) = and the projected state 
|$o) — P\^) vanishes as well. The condition that at least 
one system length must be odd is a necessary condition 
for this model. 



APPENDIX B: 

We briefly recall the transfer matrix calculation of the 
occupation probability for certain sites of the Hubbard 
chain. This calculation requires open boundary condi- 
tions. However, the influence of the boundary conditions 
decreases with increasing system size. 

The ground state of this chain can be decomposed into 
a sum over all possible configurations of N singlet bonds 
(dimers) . Every component of the ground state has ex- 
actly one dimer in each cell. Let the singlets (dimers) be 



written as b\p\ff) 



)|o). 



For every cell there exist six possible contributions: 



|1)=6L|0) |2)=6t3|0> |3)=6l4|0> 
|4)=4|0) |5)=43|0) |6)=5t3|0) 



(Bl) 



The interaction excludes any double occupancy. There- 
fore the indices a and (3 have to be different within a 
cell. This condition is not sufficient, because double oc- 
cupancy on the backbone sites could occur due to certain 
configurations of neighbor cells: If e. g. cell i contains the 
singlet |1) and cell i+1 the singlets |3), |4), or |5) the con- 
nection site would be occupied by two particles. These 
neighbor configurations have to be excluded. The restric- 
tion to allowed neighbor configurations is coded by the 
transfer matrix: 



T = 



/ 1 1 1 \ 
1 1 1 
1 1 1 
111111 
111111 

V 1 1 1 1 1 1 y 



(B2) 
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With these definitions, the norm of the ground state 
($0 1*5*0) for a chain with N cells may be written as 

(a>o|$o)=^(T^-^).,. (B3) 

The occupation of one a — 2-site well within the chain 
(the cell index m should be w 1/2N) is given by the 
expression 

(B4) 

which may be simplified by using a spectral representa- 
tion of the matrix T. The eigenvalues of T are e+ — 
^(5 + \/T3), £- = 5(5 — \/T3) and a fourfold degenerate 
eigenvalue e = 0. The right and left eigenstates may be 
calculated by r|r±) = e±|r±) and (/±|T = e±{l±\. They 
are mutually orthogonal and have to be normalized ac- 
cording to {l±\r±) = 1. Then the matrix reads 

= |r+>(6+)^(/+| + \r^){e^rm (B5) 

As usual, the contribution of the largest eigenvalue dom- 
inates the results with growing system size. In the ther- 
modynamic limit the total occupation of the a = 2-sites 
becomes 



(K,2,T 




(B6) 

The backbone occupation ma^ be calculated analogously 
or simply by applying J^a J2a=i{^z.a,a) = 2- 
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FIG. 1. Two examples of lattice structures discussed here. 
The left column shows the cells which are connected to a 
lattice (chain). Filled circles denote sites with infinite on-site 
repulsion, solid lines illustrate the hopping, a) Lattice of the 
hypercubic model for three dimensions. The sites within a 
unit cell are labeled by an index a — 1, . . . , 2D. The lattice 
consists of connected octahedra located at the vertices of a 
hypercubic lattice, b) One-dimensional Hubbard chain with 
four sites a = 1, . . . , 4 per cell. 

FIG. 2. Band structure of the non-interacting chain. The 
spectrum shows two dispersive bands separated by a gap of 
one energy unit. The third band collapses to a fiat band 
leading to a 5- peak in the corresponding density of states. 



TABLE I. Numerically calculated ground-state energies for 
systems with A'^ unit cells and 2A'' — 1 particles. The states 
are unique and correspond to fc = tt. 



N 


ground-state energy 


2 


1.2275046 


4 


1.2137261 


6 


1.2099690 


8 


1.2095401 
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